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Abstract. I describe electromagnetic model of gamma ray bursts and contrast 
its main properties and predictions with hydrodynamic fireball model and its 
magnetolrydrodynamical extension. The electromagnetic model assumes that 
rotational energy of a relativistic, stellar-mass central source (black-hole-accretion 
disk system or fast rotating neutron star) is converted into magnetic energy through 
unipolar dynamo mechanism, propagated to large distances in a form of relativistic, 
subsonic, Poynting flux-dominated wind and is dissipated directly into emitting 
particles through current-driven instabilities. Thus, there is no conversion back 
and forth between internal and bulk energies as in the case of fireball model. 
Collimating effects of magnetic hoop stresses lead to strongly non-spherical expansion 
and formation of jets. Long and short GRBs may develop in a qualitatively similar 
way, except that in case of long burst ejecta expansion has a relatively short, non- 
relativistic, strongly dissipative stage inside the star. Electromagnetic and fireball 
models (as well as strongly and weakly magnetized fireballs) lead to different early 
afterglow dynamics, before deceleration time. Finally, I discuss the models in view of 
latest observational data in the Swift era. 


THE ELECTROMAGNETIC MODEL OF GAMMA RAY BURSTS 


2 


1. Short introduction 


Gamma Ray Bursts (GRBs) are conventionally divided into two classes, short-hard and 
long-soft, distinguishe d by t heir duration (with a division near ~ 2 sec) and spectrum 
hardness (jKouveliotou et a l. 1993). Detection of Type Ic supernovae near ly coincidence 


with long GRBs unambiguously linked them with deaths of massive stars (jStanek et al 


2993tlHiorth et all 1200311 . Studies of the host galaxies of long G RBs, which turned out t o 


be actively star-forming, further strengthens this association (Diorgovski et al ., 1998) 


Recent progress in observations of short bursts showed that on one hand they show 
qualitatively similar afterglow behavior (but without any supernovae signature) while 
on the other hand their energetics was two to four orders of magnitude smaller and they 


stellar population ( 

Gehrels et al. 2995; Prochaska et a 

LO 

o 

o 

CN 

Villasenor et al ., 2995; 

Covino et al.. 2995: Rutter et al. 

, 2995; Fox et al.. 

2995 

). These indirect evidences are 


consistent with formation of short GRBs in compact star mergers (double neutron stars 


or black holes-neutron star binaries) and formation of a black hole ( e.g. iRosswog et al. 


29931: lAlov et all |2995) 


2. Short and long GRBs 

Association of two types of GRBs with different astronomical objects is somewhat 
surprising given their apparent similarity (perhaps less surprising in view of the fact that 
some short GRBs may be associated with nearby SGRs). One possible reason is that 
though short and long GRBs occur in different astrophysical setting, their appearance is 
governed by similar physical process related to formation and early evolution of stellar 
mass relativistic compact object. [Similarities o f temporal and spectra l properties of the 


2992; 


Ghirlanda et al. 


first 2 seconds in long bursts and short bursts (INakar & P iran, 

2993) may be an indication of this.] But then one expects that during merger of, e.g. 


, two neutron stars the resulting black hole has large angular momentum and thus can 
potentially release much more energy than observed (one can invoke different efficiencies, 
but a naive guess would be that it’s harder to extract and propagate energy from a 
compact object inside a stellar core, contrary to observations). 

So why short GRBs are so under-energetic if compared with long ones? One 
possibility is that presence of a disk is a necessary condition for extracting energy from 
a black hole, so after disk disappears energy extraction stops. In this case the energy 
and angular momentum that will power a GRB outflow effectively come not from the 
central black hole but from the surrounding disk (in a sense that its the energy and/or 
lifetime of the disk and not the energy in the black hole tha t det ermines the resulting 


energy of a GRB outflow) (for related discussion see also v an Putten , 2995) 


In case of neutron star mergers, black hole forms fairly early, while the mass of the 
accretion disk is small, < 9.1 M 0 , with short viscous time scales, ~ 9.1 — 1 sec (e.g. 
Ruffert & Janka, 2991). On the other hand, a black hole inside a collapsing core of 



















































































THE ELECTROMAGNETIC MODEL OF GAMMA RAY BURSTS 


3 


a massive star may accrete several solar masses of material {e.g. iMacFadven et al. 
2001) (at any given moment the mass of the disk is small, but large amount of 
mass, ~ 1 — 10 M q passes through the disk during accretion). In addition, amount 
of the rotational energy stored in case of core collapse depends on core rotation before 
explosion (which, in turn, depends on metallicity through wind angular momentum loss 
Wooslev & Heger ( 20050 ). resulting in broad spread of rotational energies. 


3. Principal issues: electromagnetic and fireball models 

In this contribution we describe electromagnetic model of GRBs which assumes that the 
energy that will power a GRB comes from rotational kinetic energy of a central source. 
The energy is extracted t hrough magnetic field, which can be generated by a local 


dynamo mechanisms (e.g. iDe Villiers et all 120051: lHawlev fc Krolikl. l2005t iProga et oL. 
2003 : McKinney fc Ga/mmiel. 2004 h As argued above, the GRB energy should then be 


related not to the total rotational energy of a central black hole but to the disk around 
it. Another possibility for long bursts is formation of a ”millisecond magnetar”, a fast 
rotating strongly magnetized protoneutron star. 

Whether magnetic fields play an important dynamical role at any stage in the 
outflow remains, in our view, one of the principal issues in GRBs physics. Currently, the 
overwhelming point of view, advocated by the fireball model (FBM) is that magnetic 
fields do not play any major dynamical role (except, perhaps, at a very early stage; 
after which fields are dissipated quickly). FBM advocates that in the emission region 
m agnet i c fie l d are re -crea ted locally {e.g. through development of Weibel instability 
( Medvedev & Loebl 1999 )). with energy density typically much smaller than plasma 
energy density. Fields are small scale, with correlation length l c much smaller than 
the ’’horizon” length l c -C R/T {R is the radius of the outflow in the laboratory 
frame and F is its bulk Lorentz factor). Alternative approach, advocated by M HD 
and e lectromagnetic models {e.g. lUsovL 11992c iBlandfordl 120021: iLvutikov & Blandford , 
2003) is that there are dynamically important large scale fields with ’’super-horizon” 
correlation length l c > R/T , which are created at the source and which may play a 
major role in driving the whole outflow in the first place. 

To quantify the dynamical importance of large scale magnetic fields, it is useful to 
introduce magnetization parameter a as a ratio of Poynting Fp oynt i ng to (cold) particle 
F p fluxes (or as a ratio of rest frame energy densities) 

Apoynting B 2 b' 2 


a = 


7 p 47rr pc 2 87ip'c 2 


( 1 ) 


where B and p are magnetic field and plasma density in the lab frame, b' and p' are 
magnetic field and plasma density in the plasma frame (where electric field is zero). 
For (T < 1 magnetic fields are dynamically unimportant (this is assumed within a 
framework of a conventional FBM), while for a > 1 magnetic fields start to play 
important dynamical role. For cr 3> 1, there is an important qualitative change in 
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the dynamical behavior of the flow at a cr n = T 2 /2. For cr < a crit the flow is super- 
Alfvenic, while for a > a crit the flow is sub-Alfvenic (and sub-fastmagnetosonic). The 
different is somewhat analogous to the difference between sub-sonic and supersonic flows 
in hydrodynamics. 

Thus, depending on the parameter a three qualitatively different regimes for 
expansion of the ejecta may be identified, which I will call (i) fireball model (FBM 
below) a -C 1, (ii) MHD models 1 < a < a cr i t and (iii) electromagnetic model (EMM 
below) cr > a cr it . These three possibilities leads to a qualitatively different dynamic 
behavior of flows. Let us next describe qualitatively main features of the models. (As 
the FBM and EMM are at the extre me rang e of cr we discuss those first.) f 


Fireball model (FBM, e.g. 


Piran 


2004): The defining characteristic of the 


FBM is that at intermediate distances (far from the central source but before most 
energy is transfered to the forward shock) most energy produced by the central source 
is carried by the bulk motion of ions. In temporal order the transformations of energy are 
as follows. Initially, the energy that will power the GRB and its afterglow is thermalized 
near the central source, so that most of it is converted into lepto-photonic plasma. This 
internal energy is then converted to the bulk motion of ions, and reconverted back into 
internal at internal shocks; at the same time, small scale magnetic fields are generated. 
The energy of these generated magnetic fields is then used to accelerate leptons via 

Fermi mechanism to highly relativistic ener gies §. _ 

Electro-magnetic model (EMM, Lvutikov &: Blandford . 2003 1. The 
defining characteristic of the electro-magnetic model is that the bulk energy of the 
flow is carried subsonically by magnetic field. In temporal order the evolution of the 
energy proceeds as follows. The energy that will power a GRB comes from kinetic 
rotational energy of the central source (millisecond pulsar or BH-disk system). It is then 
converted to magnetic energy using unipolar inductor (like in pulsars), transported to 
large distances in a form of strongly magnetized wind and is used to accelerate particle 
in the emission region. Acceleration of particles is done via magnetic dissipation (not 


MHD model (e.g. Drenkhahn Spruit, 

2002 

): In this case most energy 

is also curried by magnetic held (similar to Lvutikov & Blandford (20031)), but the 

how is supersonic, similar to FBM. In the current version ( 
magnetic held energy is hrst converted into bulk motion 
internal shocks, similar to FBM. This is, in principle, not 

Drenkhahn & Spruit. 2002), 
and then dissipated through 
necessary, so that magnetic 


held energy may be dissipated directly into emitting particles, similar to EMM. 

The principal differences between EMM and MHD approaches is that MHD-type 
outflow usually cross fast magnetosonic cr itical surface after which moment they become 
causally disconnected from their source (Goldreich & Julian, [1970). Initially the how 


$ Definitions and discussion below is based not on the nature of the central object but on ejecta content 
at large distances from the central region and before production of y-rays occurs. 

§ The energy that goes to non-thermal particles is electro-magnetic even in the FBM: Fermi-type 
acceleration is done by turbulent EMF associated with fluctuations of magnetic field 
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is expanding freely, so that the flow dynamics is determined by the internal structure 
of the flow. Only after the flow reaches the terminal velocity, the interaction with the 
mediums becomes important. Unlike MHD, force-free flows are sub-fastmagnetosonic 
so no conditions at the fast critical surface appear. In this case it is the interaction with 
boundaries that determines the properties of the flow (similar to subsonic hydrodynamic 
flows). Thus, the distinctive feature between MHD and force-free flows is whether 
the wind becomes fast supersonic (MHD regime, a < cr crit ) or not (force-free regime, 
cr > cr crit ). This important difference leads to somewhat different dynamics of the flow 
and can be tested with observations, as discussed in 


4. Source Formation and Energy Release in EMM 


4-1. Electromagnetic luminosity and currents 

In this section I will describe main ingredients of the EMM, stressing its principal 
difference and predictions from the FBM. EMM assumes that the GRB “prime mover” 
is relativistic, fast rotating, near stellar-mass source. As discussed above, in order 
to reconcile energies of short and long GRBs the “prime mover” should be not the 
black hole but the disk around it. For numerical estimates I will assume that a central 
source generates luminosity L = L 50 x 10 50 erg/s for a time t s , where t s ~ 100 s 
for long bursts and t s ~ 1 s for short ones (there are indications that both long and 


short bursts are powered by the same luminosity, but for different time (Fox et al 


2005)). The mass of the central source is ~ 0.01 M & for short bursts and 


Mr 


0 


for long bursts (we stress again that this is the total mass passing through the disk, 
not in the black hole). The source is assumed to rotate with a spin frequency 
~ kHz. In addition, it is assumed that a source possesses a large magnetic field of 


Wooslev & Hegerl 2005), 


B s ~ 10 14 L^ 2 G. If initially the core is fast rotating ( e.g. 
the total rotational energy in the disk, ~ Mdi S kR d isk ^ 2 , i n case °f core collapse is much 
larger, ~ 7.9 x 10 52 {Mdi S k/M Q ){Rdi S k/ x 10 6 cm) 2 (D/6.28 x 10 3 rad/s) 2 erg, than in the 
case of mergers, ~ 7.9 x 10 5 ° erg with M disk ~ 0.01 M q in the above estimate. The 
source is expected to be active for t s ~ E/L ~ 100 s for longs and ~ 1 s for shorts. 

Rotational energy of the central object is extracted by magnetic fields through 


unipolar induction mechanism, similar to prevailing model of AGN jets {e.g. iFerrari 
120041) . Magnetic fields both launches the jet {e.g. through Blandford-Znajek mechanism) 
and collimates it by hoop stresses. Latest full relativistic MHD numerical simulations of 
accretion dis k-black hole systems do show formation of the strong ly-magnetized axial 


funnel {e.g. De Villiers et al 


2005 


McKinney & Ganunie, 2004). Thus, large scale, 


energetically dominant magnetic fields may be expected in the launching region of GRB 
jets. 

Qualitatively, in the immediate vicinity of the source, the plasma is separated into 
two phases: an internal matter-dominated phase in which large currents are flowing 
and external magnetically dominated phase. Strong magnetic fields and magnetic flux 
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are generated in the dense medium (a disk or a differentially rotating neutron star-like 
object), while relativistic outflow is generated in the magnetically dominated phase. In 
this case matter loading may be expected to be small (e.g. analogous to pulsar wind). 
As the source remains active for ~ thousand to million dynamical times the flow will be 
able to settle down quickly to a quasi-steady state evolving slowly as the hole or neutron 
star slows down. The separation into matter- and magnetic field-dominated phases is 
somewhat similar to the Sun, where dynamo operates in tachocline, deep below the 
surface, while magnetic energy, and most importantly magnetic flux, are dissipated 
outside the star. 

The key assumption of the model is that the dissipation rate at the source remains 
low enough so that the power continues to be dominated by the electromagnetic 
component (rather than the heat of a fireball) well out into the emission region. Thus, 
electrical currents flows all the way out to the expanding blast wave, rather than been 
dissipated close to the source. For electromagnetically dominated outflows the value of 
the total current may be related to the total luminosity of the source 

1 ~ ~ 3 x 102 ° L 50 2 A > ( 2 ) 

where notation X n = (Ai/10 n ) was adopted. Under general electromagnetic and 
relativistic conditions the total impedance of the source and the emission region is 
close to the impedance of free space Z ~ 100 0 

In case of long burst, associated with collapse of massive stars, the source will 
initially inflate a non-relativistically expanding electromagnetic bubble inside the star. 
This magnetized cavity is separated from the outside material by the (tangential) 
contact discontinuity (CD) containing a surface Chapman-Ferraro current. This current 
terminates the magnetic field and completes the circuit that is driven by the source. On 
a microphysical level the current is created by the particle of the surrounding medium 
completing half a turn in the magnetic field of the bubble, so that the thickness of the 
current-currying layer is of the order of ion gyro-radius. After the break out the density 
that controls the ejecta expansion falls down, so that expansion becomes relativistic. In 
case of short burst, associated with merger of two neutron stars, a somewhat similar 
process will happen, except that there is no non-relativistic stage, so that the bubble is 
directly inflated in the circumburst medium. 


4-2. Distribution of current in the wind: structured jet 


The form of expanding bubble depends on lateral distribution of source luminosity, 
which within a framework of EMM mode is given by lateral distribution of current. At 
relativistic stage expansion is nearly ballistic (Shapiro, 1979) while at a non-relativistic 
stage inside a star a flow may be collimated both through the action of magnetic hoop 
stresses and interaction with the surrounding gas (see H5.2j) . One particular stationary 
outflow configuration which captures the essential features of the outflow, is that the 
outgoing current is confined to the poles and the equatorial plane and closes along the 
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surface of the bubble, Fig. IA11 This current distribution minimizes the total energy 
given a total _toroidal m agnet ic flux and has been advocated in relativistic stationary 
winds ( Hevvaerts Hr. Norman! 2003). The magnetic held in the bubble is inversely 
proportional to the cylindrical radius, B^ oc l/^rsin#). Accompanying this magnetic 
held is a poloidal electrical held so that there is a near radial Poynting hux carrying 
energy away from the source. Thus, as the magnetic held strength is strongest close 
to the symmetry axis, the bubble will expand fastest along the polar direction. The 
i nternal structure of a n outhow then cor responds to a ’’structured jet” with Lg ~ 0 2 


flLi punov et al. 


2001 


Rossi et all .2002), so that the central source releases an equal 


amount of energy per decade of 6. 


5. Long GRBs: expansion inside a star 

5.1. How important is dissipation? 

Some fraction of the central source luminosity is likely to be dissipated close to the 
source. The hreball model implicitly assumed that all of the energy released is quickly 
converted into heat, while in the electromagnetic model this does not happen and the 
energy hows way from the light cylinder mainly in the form of an electromagnetic 
Poynting hux and the load impedance is located in the emission region. 

The issue of dissipation is somewhat complicated as I discuss next. On one hand, 
somewhat paradoxically, it becomes harder to convert electromagnetic energy directly to 
pair plasma the stronger the magnetic held becomes. The reasons is that the maximum 
potential drop that is available for dissipation will be limited by various mechanisms of 
pair production. Typically, after an electron has passed through a potential difference 
AV ~ 10 9 — 10 12 V it will produce an electron-positron pair either through the emission 
of curvature photon or via inverse Compton scattering. This will be followed by an 
electromagnetic cascade and the newly born pairs will create a charge density that 
would shut-off the accelerating electric held. Put another way, the pair density required 
to supply the electrical current and space charge scales linearly with the held strength, 
while the electromagnetic energy density scales as its square. The stronger the held, the 
more likely it is to persist into the outhow. It is because GRBs are so powerful that the 
dissipation in the source is probably low. 

There is an important caveat to the preceding discussion which applies to very early, 
nonrelativistic, stages of bubbl e expans ion in case of long bursts. Somewhat similar 
to pulsar wind nebulae (c/. [Rees fc Gunnl . Il974 h non-relativistic, ideal, homologous 
expansion of strongly magnetized nebular which is causally disconnected from the source 
and injects toroidal magnetic field with nearly a speed of light cannot occur. The reason 
is that magnetic hux and energy are supplied to the inflating bubble by a rate that 
cannot be accommodated in the bubble. The rate of supply is determined by the 
processes inside the light cylinder of a newly-formed, compact object, while inflation of 
the bubble is controlled by the external gas density. (Even if the wind remains subsonic, 















THE ELECTROMAGNETIC MODEL OF GAMMA RAY BURSTS 


it is unlikely that processes at the edge of the inflating bubble would influence wind 
generation region near the light cylinder.) This leads to the following ’’contradiction”. 
Magnetic flux (integrated over the meridional plane) is supplied to the bubble at a 
rate ~ 21 c. Similarly, energy is supplied to the bubble at a rate Uem ~ SI. We 
can also compute the magnetic flux <f> = Cl and the energy stored within the bubble 
Uem = SI 2 1 2, using the self inductance C. Let the bubble radius be R{9R) where 9 is 
the polar angle measured from the symmetry axis defined by the spin of the compact 
object. If the magnetic field in the bubble is predominantly toroidal between cylindrical 
= R sin 6 , this is given by 

[dzln 

2tt J \ uj _ 

We therefore see that, if the bubble expands sub-relativistically, the rate of supply of 
both flux and energy exceeds the rate at which the flux a nd ene rgy can be stored by 


radii zu rri j n and zu max 


C 


( 3 ) 


a factor ~ [In (zj mSLX /w ia i n )(dz/dt)/c] 1 (c/. Rees & Gunn], 19741). Therefore, for non- 


relativistic expansion of the bubble too much flux and too much energy is generated by 
the source. 

The way out of the ’’paradox” is that dissipation must become important that 
will destroy some magnetic energy and most importantly eliminate most of the toroidal 
flux. Most of the dissipation is likely to occur near the axis where the current density is 
highest and the susceptibility to current-driven instability is the greatest. In this case a 
lateral flow of energy will set in carrying the poloidal field lines with it towards the axis. 
This, in turn, will leads to the pilc-up of magnetic field near the axis and to f aster radial 
expansion near the axis (the toothpaste tube effect) ( Lvutikov fc Blandfordl. 2003). 


5.2. Form of the expanding bubble 

The dynamics of a non-spherically expanding bubble may be described using the method 
of Komnaneetsl (J1960). Consider a small section of non-spherical non-relativistically 
expanding CD with radius R(t,9). The CD expands under the pressure of magnetic 
held so that the normal magnetic stress at the bubble surface is balanced by the ram 
pressure of the surrounding medium. At the spherical polar angle 9 the CD propagates 
at an angle tana = — dl g g R to the radius vector. Balancing the pressure inside the 
bubble B 2 /(8tt ) = I 2 / (2ttc 2 R 2 ) with the pressure of the shocked plasma gives 



= k- 


m 


i + 


'5 In R N 
89 , 


( 4 ) 


2t:R 2 sin 2 9p(R, 9) 

where k is a coefficient of the order of unity which relates the pressure at the CD to the 
pressure at the forward shock. 

Equation (ED shows that non-spherical expansion inside the star is due both to the 
anisotropic driving by magnetic fields and collimating effects of the stellar material (the 
term in parenthesis, which under certain conditions tends to amplify non-sphericity). 
The rate of expansion of the bubble inside the star depends upon the density profile of 
the stellar envelope and the time evolution of the luminosity (or, equivalently, of the 
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current /(f)). For a given dependence p(R,9) and I(t) Eq. (HJ) determines the velocity 
of the CD. Generally solutions will be strongly elongated along the axis. A simple 
analytical solution for /, p ~ const is 


RM = (~— 2 

\7T pC Z 


1/4 


Vt_ 

sin 9 


(5) 


(current is related to the luminosity by Eq. (j2J). Qualitatively, the bubble and the 
forward shock will cross the iron core (r c ~ 2.5 x 1 0 8 crril in t ~ rywW B (rA ~ .3 sec, 
short enough to produce an ample supply of 56 Ni (jVVooslev et al . 2003 ). If M(R ) is 
the stellar mass external to radius R, then the breakout time is 


tbreakout (0) 


ie\ 


M 


— , V2 (^X /2 l~^ 
M q ) UoJ * 


7 50 


( 6 ) 


The electromagnetic bubble can be confined equatorially by the star for the duration of 
the burst f breakout^/ 2 ) 100s and will expand non-relativistically as we have assumed. 
However the expansion along the axis proceeds on a short time scale and breakout should 
occur early in the burst. 

Thus, along the jet axis non-relativistic expansion lasts for several seconds, which 
is much shorter than the burst duration. After breakout the flow quickly becomes 
relativistic so there is no need for dissipation anymore. Thus, relative fraction of 
dissipated energy is small, so that overall the flow magnetization remains large. Most of 
the dissipation described above will result in creation of lepto-photonic plasma, which 
decouples after photosphere, so that the remaining flow remains strongly magnetized. 
In summary, when the bubble expands non-relativistically, it must be dissipative, while 
after breakout, the expansion becomes relativistic and the resistance falls so that the 
electromagnetic energy that is still being supplied by the source is mostly absorbed by 
the inflating bubble and by doing work against the surroundings. 


6. Optically thick expansion: mini-fireball 

Under the electromagnetic hypothesis, most of the energy released by the source comes 
out in the form of Poynting flux. However there must be some dissipation that would 
lead to creation of a lepto-photonic component (as discussed in m, in case of non- 
relativistic stage of bubble expansion inside a star dissipation may be considerable). 
This will create an optically thick ’’warm” fireball (in a sense that it is dominated by 
magnetic field energy, but also has considerable pressure). Expansion of this ’’warm” 
fireball will create a thermal precursor, similar to the conventional FBM, but modi fied 
by presence of magnetic field (ILvutikov fc Uso\i. 12000c iLvutikov fc Blandford . 2003). 

At early stages (before the breakout in case of long bursts) the plasma enthalpy is 
strongly dominated by lepto-photonic plasma with a temperature 


T 


( _ h _ 

l aADr 2 /3cT 2 (l + cr) 


1/4 


(7) 
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where All is a typical opening solid angle and luminosity of the source then can be 
written as L = / dLlT 2 r 2 /3c (b 2 /2 + w), where w is plasma enthalpy, b is a toroidal 
magnetic held in the plasma rest frame times \/4vr, 

After breakout the how will accelerate to relativistic velocities. Initially, conical 
expansion is mostly pressure-driven, even in the strongly magnetized case (in this case 
magnetic pressure gradients and hoop stresses balance out each other). This results in 
dynamics qualitatively similar to the unmagnetized case: the wind plasma accelerates 
T ~ r while its density, pressure and temperature decrease n ~ r~ 3 , v ~ r~ 4 . T ~ r" 1 . so 


that magnetization parameter remains approximately constant (Lvutikov & Blandford 

<H)). 

When the temperature falls below ~ 10 —20A;eU, most of the pairs annihilate. This 
suddenly reduces the optical depth to Thomson scattering below unity. (Under certain 
conditions photons may remain trapped in the how. In this case, thermal driving by 
photon pressure continues, until the thermal photons escape.) As a result the lepto- 
photonic part of the how decouples from the magnetic held and cr increases by roughly 
seven orders of magnitude to a ~ 10 9 . The thermal radiation from the lepto-photonic 
component has a rest-frame temperature T 0 ~ 10 — 20 keV times a boost due to the 
bulk motion. This thermal radiation, which s hould peak around ~ 100 keV may put 


constraints on the initial cr (Lvutikov & Usov, 


2000 ; 


Daigne & Mochkovitch. *2002 


7. Relativistic expansion 

7.1. Short GRBs and long GRBs after breakout 

In case of long GRBs inflating a bubble inside a star, eventually the bubble will break 
free out of the star forming two axial jets along which Poynting hux will how until the 
central source slows down on the time scale t s ~ 100 s. Outside the star, the bubble will 
expand ultra-relativistically and bi-conically. For short GRBs, presumably associated 
with merger of neutron stars in a low density environment, there is no preceding non- 
relativistic stage so expansion of the bubble is relativistic from the beginning. In case of 
relativistic motion there is no necessity anymore to destroy magnetic hux through ohmic 
dissipation: the effective load can consist of the performance of work on the expanding 
blast wave. This is where most of the power that is generated by the central magnetic 
rotator ends up. 

After the bubble has expanded beyond a radius r sh rs./ ct s ~ 3 x 10 12 G ]2 cni 
(~ 10 10 cm for short bursts) the electromagnetic energy will be concentrated within 
an expanding, electromagnetic shell with thickness ~ r s h and with most of the return 
current completing along its trailing surface (see Fig. EH). The global dynamics of this 
shell and its subsequent expansion are set in place by the electromagnetic conditions at 
the light cylinder and within the collimation region. An im por tant pr operty o f ultra- 


relativistic outflows is that they are hard to collimate (jGhiueh et a l. 1991; iBogova lov 


2001), so that any collimation should be achieved close to the source, within a star, 
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where the flow is only mildly relativistic. 

Interaction of the magnetic shell with the circumstcllar medium proceeds in a 
similar way to the non-relativistic expansion inside a star: the leading surface of 
the shell is separated by a contact discontinuity (which actually becomes a rotational 
discontinuity if the circumstellar medium is magnetized (Lvutikov, 2002a)). Outside 
the CD an ultra-relativistic shock front forms and propagate into the surrounding 
circumstcllar medium. The expansion is non-spherical. As long as the outfl ow is ultra- 
relativistic, the motion of the forward shock is virtually ballistic (Shapiro, 111 970 ) and 
determined by the balance between the magnetic stress at the CD and the ram pressure 
of the circumstellar medium. 

A type of collimation in case of electromagnetic explosions is somewhat different 
from the conventional jet models of AGNs. We expect that large Poynting fluxes 
associated with explosive release of ~ 10 51 ergs in case of long GRBs are sufficient 
to drive a relativistic outflow over a large solid angle, so that during the relativistic 
stage the resulting cavity is almost spherical, but the Lorentz factor T of the CD is 
a strong function of the polar angle. The angular distribution of magnetic field (and 
of the Lorentz factor of the expansion) depends on the dynamics of the bubble at the 
non-relativistic stage and the distribution of the source luminosity. 

In the framework of electromagnetic model the outflow is strongly magnetized and 
subsonic, a > a cr u , despite been strongly relativistic. In this case the ejecta, in some 
sense, may be considered as a collection of outgoing fast magnetosonic waves propagating 
from the source to the contact discontinuity. Motion of the CD is then determined by 
the pressure balance between the Poynting flux from the source and the ram pressure of 
the ISM. Thus, motion of the CD depends on the source luminosity L(t') at the retarded 
time t' such that R(t ) = t — t'. In addition to forward flux, there is a much weaker, by a 
factor T 2 , reflected flux that propagates backward into the flow the information about the 
circumstellar medium. Interference of forward and backward propagating waves allows 
to define a finite Lorentz factor of the ejecta. The distribution of reflected current 
is determined by the outgoing current and the boundary conditions. At later times 
multiple reflections from the contact discontinuity and the center become important as 
well. 


1.2. Stages of relativistic expansion of electromagnetic shell 

Relativistic expansion of the magnetized shell may be separated into two stages, which 
we will call ’’early” and ’Tate”, depending on whether or not most of the fast waves 
emitted by the central source have caught up with the CD and their energy has been 
given to the circumburst medium. The transition between two stages occurs at the 
moment which is similar to the deceleration radius in fireball model, except that in the 
case of EMM the shell is decelerating all the time, but with different laws before and 
after the transition. Keeping with the tradition we will still call the transition radius as 
deceleration radius. 
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7.2.1. “Early” stage. At r > r ph the outflow becomes a relativistically expanding 
shell of thickness ~ cf s ~ 3 x 10 12 cm for long GRBs and ~ 3 x 10 10 cm for short 
GRBs. The shell contains toroidal magnetic field; the current now detaches from the 
source and completes along the shell’s inner surface. At this stage the CD is constantly 
re-energized by the fast-magnetosonic waves propagating from the central source. The 
average motion of the CD R(t) is determined by the average luminosity at the retarded 
time t': 


L n (t') ~ pc 3 T 4 R(t) 2 f3 3 ( 8 ) 

which for constant luminosity gives T rsj (■ Ln/pc 3 )V 4 r-V 2 (in a constant density medium) 
or T ~ (Lq./AkpqTqC 3 ) 1 / 4 = const (in a p(r) = Po( r o/' r ) 2 wind). If the central source 
releases most of the current along the axis and the equatorial plane, as argued in £ 14.21 
then T oc 1 / \J sin 0 at this stage. If source’s luminosity varies, this will be reflected in 
the ’’jitter” of the CD. Development jT i nstab ilitie s at the CD, like impulsive Kruskal- 
Schwarzschild instability ({L yu tikov & Blandford . 12003 ) may lead to dissipation and 
particle acceleration. The internal structure of the magnetic shell is a messy mixture of 
the outgoing waves from the source and the ingoing waves reflected from the CD, similar 
to a pre-Sedov phase in hydrodynamical explosions. Unlike the case of a hydrodynamic 
blast wave with energy supply, no internal discontinuities form inside magnetic shell. 

Early stage lasts for ct s < r < r dec , where 

r dec = ( L n t 2 Jpc) 1/4 ~3x 10 16 L^ 4 t^ 2 n _1/4 cm, (9) 


for long bursts (in the observer frame this phase lasts ~ t s ~ 100 s) and r dec ~ 3 x 10 15 cm 
for shorts (similarly, in the observer frame this phase lasts ~ t s ~ 1 s). Radius r dec (JUJ) 
is somewhat similar to the deceleration radius in case of FBM; at this moment most 
energy of the shell in given to the circumburst medium, Lot s rv./ E(0) ~ pc 2 r 3 dec T(r dec ,6) 2 
(note that here T = T(r dec ), not T 0 as in case of FBM, since there is no formal definition 
of Tq in case of EMM). 


7.2.2. ’’Late stage” (r dec < r < tnr = (L n t s / pc 2 ) 1 ^ 3 ). At distances r > r dec most of the 
waves reflected from the CD have propagated throughout the shell, so that all the regions 
of the shell come into causal contact. Most of the energy of the explosion will reside in 
the blast wave which will eventually settle down to follow a self-similar expansion. As the 
expanding shell performs work on the surrounding medium its total energy decreases; 
the amount of energy that remains in the ejecta shell during the late stage is small, 
~ Efi/T 2 . Most of the energy is still concentrated in a thin shell with A R ~ R/T 2 near 
the surface of the shell which is moving according to T ~ Eq/ pc 2 r~ 3 ^ 2 (in a constant 
density medium), or T ~ r^ 1 / 2 (in ap~ r~ 2 wind). If the central source releases most of 
the current along the axis and the equatorial plane, as discussed in £14.21 then T oc 1/ sin 6 
at this stage. [Note that at the ’’early stage” T oc l/vsind, but no lateral re-distribution 
of energy is required at the transition since transition between ’’early” and ’’late” stages 
occur at different times for different 9] The energy of the shell decreases only weakly 
with radius, dE/ctt ~ 1/r in constant density and dE/dt ~ 1/r 3 in a wind, so that the 
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surface of the shell keeps moving relativistically as long as the preceding shock wave 
is moving relativistically, until r ~ (Ep/pc 2 ) 1 ^ 3 ~ 10 18 cm - the shock never becomes 


completely free of the shell ( Lvutikov & Blandfordl . 2003). Interestingly, the structure 
of the magnetic shell (in particular the distribution of energy) resemble s at t his stag e 
the structure of the hydrodynamical relativistic blast wave (Blandford & McKee, 1976). 
This can be formally understood by noting that for motion perpendicular to magnetic 
field dynamical equations for magnetized f low can b e reduce to non-magnetic case, but 


with a different equation of state (jLandau fc Lifshit.d . 1 19821) . 


” Late stage” of magnetic shell expansion corresponds to the conventional afterglow 
phase when synchrotron and inverse Compton radiation is emitted throughout the 
electromagnetic spectrum. The initially aspheric expansion will give the appearance of 
a jet with the “achromatic break” occurring when the fastest Lorentz factor of the spine 
T{0 = 0 ) becomes comparable with the reciprocal of the observer’s inclination angle 
with respect to the symmetry axis, T (0 = 0 ) rs./ l/9 ob . When r > r NR r s-' (Lf/pc 2 ) 1 / 3 ~ 
2 x 10 1S lIq 3 ^ 3 n^ 1//3 cm, the blast wave become non-relativistic and will become more 
spherically symmetric, while evolving towards a Sedov solution. 


8. Production of GRB 

By the time the shell radius expands to r dec most of the electromagnetic Poynting flux 
from the source will have caught up with the CD and been reflected by it, transferring 
its momentum to the blast wave. Simultaneousl y a str o ng regi on of magnetic shear is 
likely to develop at the outer part of the CD (ILvut i kovi 2002a). 

We propose that the 7 -ray-emitting electrons are accelerated near r dec ~ 10 15 — 
10 16 cm (for long bursts) and ~ 10 15 cm (for short bursts) due to development of 
electromagnetic current-driven instabilities (conventional model of particle acceleration 
- acceleration at internal shocks - cannot work in this model since in the limit a 1 fast 
shocks are either weak or do not form at all). The development of current instabilities 
usually results in enhanced or anomalous plasma resistivity which leads to an efficient 
dissipation of the magnetic held. The magnetic energy is converted into heat, plasma 
bulk motion and, most importantly, into high energy particles, which, in turn, are 
responsible for the production of the prompt 7 -ray emission. The conversion of magnetic 
energy into particles may be very efficient. For example, recent RHESSI observations 
of the Sun indicate that, in reconnection regions, most magnetic energy goes into non- 


thermal electrons (jBenz et all 120031) . 


To illustrate how magnetic dissipation may proceed, we describe shortly physics 
underlying the development of the so-called tearing mode. Consider a smooth 
distribution of electrical current, which can be viewed as a set of many small current 
wires. Since parallel currents attract, such a system is likely to develop narrow current 
sublayers where dissipation, which is inversely proportional to the square of magnetic 
held gradient, becomes high. In addition, high anomalous resistivity, proportional to 
local current density is likely to develop. Similar to non-relativistic plasmas, in strongly 
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magnetized plasmas tearing mode develops on time scales much shorter than resistive 


time scale in the bulk (Lvutikov, 


200 


3 ) 


The final outcome of the development of the 


tearing mode is formation of reconnection sites and dissipation of magnetic energy. 

Particle acceleration by dissipative magnetic fields may proceed in a number of 
ways. The best studied non-relativistic example is particle acceleration in reconnection 
regions either by inductive electric fiel ds ou tside the current sheet or resistive electric 


fields inside the current sheets ( e.g. ICraig fc Litvinenkol. 120021) or format ion of shocks in 


the downstream of reconnection regions (e.g. Blackman fc Fieldl 1004 b Investigation 
of th e par ticle acce leration in the relativistic regime of reconnection is only beginning 
Larrabee et al . 20021). Relativistic reconnect i on may produce powe r-law spectra of 


[e.g- - - - - , r - - -x-- 

accelerated particles (IZemtani and Hoshinol 12004 lLarrabee et al. . 2002 ). For exam ple, 
in the relativistic Sweet-Parker reconnection model (Lvutikov & Uzdenskv, 2003), if 
one balances linear acceleration inside the reconnection layer by the resistive electric 
field, d t £ ~ eEc with the rate of particle escape (proportional to relativistic gyro- 
frequency), d t lnN (£) ~ u>B(mc 2 /'£), one finds N(£) ~ w here £ j s the 

energy of a partic l e, N(£) is the particle number and (3i n is the inflow velocity 
(Zenifani and Hoshino, 2004)- For relativistic reconnection the inflow velocity can be 
relativistic (Lvutikov & IJzdenskvI. 1200311 . fl in —> 1. The fact that reconnection models 
can produce spectra which are prohibitively hard for shock acceleration may serve as a 
distinctive property of electromagnetic models. 

In addition to the acceleration mechanisms which are based on known non- 
relativistic schemes, it is feasible that acceleration in relativistic, strongly magnetized 
plasma may proceed through mechanisms that do not have non-relativistic or fluid 
analogues. Examples of this type of acceleration include particle acceleration through 
formation of a spectral cascade of nonlinear waves in force-free plasma which transfer 
energy to progressively larger wave vectors until this energy is take n up in accelerating 
a population of relativistic electrons and positrons Thompson & Blaes ( 19981 ). Since 
for cr > 1 cascade is likely to be terminated at plasma frequency, which is lower by a 
factor yfo the cyclotron frequency, the likely emission mechanism in this case is inverse 
Compton scattering. Another possibility is development of kinetic electroma gnetic -type 
instabilities of the shell surfac e currents, a s prop osed by ISmolskv fc IJsovI 1)19961) and 
in somewhat different form by Liang et al. (12003)1) . Since studies of kinetic properties 
of strongly magnetized relativistic plasmas are only beginning, it is hard to predict 
acceleration efficiency and particle spectra. Numerical studies in the coming years will 
be most important here. 


9. Production of afterglows 

Except at early stage (as discussed in U7.2.1D . afterglows are generated in a similar 
way both in the FBM and EMM. As the magnetic shell expands, its energy is 
gradually transfered to the preceding forward shock wave. Relativistic particles are 
accelerated in the blast wave producing the observed afterglow in a manner which is 
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similar to what is proposed for fluid models, except that contact discontinuity itself 
may be an im portant sou r ce o f ma gnetic flux through impulsive Kruskal-Schwarzschild 
instability (jbvntikov fc Bla/ndfordl 120031) . so that afterglows may result from a mixture 
of relativistic particles, derived from the shock with magnetic held derived from the 

shell. 

At late times, well beyond r^c (which in observer frame is nearly coincident with 
the prompt phase), the temporal behavior of proper afterglow (as opposed to tails of 
prompt emission, see mm is determined by the total energy release Eq and not by the 
form of that energy. As a consequence, late afterglow observations can hardly be used 
to distinguish between the models. The only property of the source that the forward 
shock ’’remembers” at late times is the angular distribution of the deposited energy 
E{6 ) (there is little sideways evolution in relativistic regime (Sha p iro. 19791 1). Thus, 
the angular distribution of the total energy E{6) can be used to distinguish between 
different models, if a model predicts it. 

In case of EMM, the preferred lateral distribution of the magnetic held, energy 
fluxes and luminosity correspond to the line current, H4.21 so that L l/sin 2 #. This 
translates into distribution of Lorentz factors of the forward shock 


T 


E 


1/2 


t-3/2 


,PexC 5 J + 01 


where 6q is the angular size of the core of the jet 

which gives 9 0 


( 10 ) 

(Its minimal size is magnetic 


ueoye racuus, rn 

Lvutikov & Blanc 

- v 

ford. 

q WJ 

2003).) 

(or universal jet) ( 

Lipunov et al. 


m? c s a 2 r 2 
Le 2 


1/4 


m_3r- i/4 i /2 r 1/2 
IU ±J 50 a 9 I 2 1 

This type of shock has b een named “structured jet” 


2001 


Rossi et all 12002 ). though in our model there 


is no proper “jet”, but simply a non-spherical outflow. The most intense bursts and 
afterglows in a flux-limited sample will be seen pole-on and can exhibit achromatic 
breaks when T ~ 6~ 1 , which might be mistaken for jets. 

In conclusion, observational appearance of GRB afterglows depends mostly on two 
parameters: (i) explosion energy (more precisely, on the ratio on the explosion energy 
to circumstellar density) and (ii) the viewing angle that the progenitor’s axis is making 
with the line of sight. This possibility, that all GRBs (and XRFs) are virtually the same 
but viewed at different angles resembles unification scheme of AGNs. 


10. Tests of GRB models 

10.1. Testing the fireball model: reverse shock emission 

Perhaps the simplest test of GRB models could have come from observations of 
emission from the reverse shock propagating in the ejecta, which typically falls into 
the optical range. FBM predicts strong reverse shock emission, so that absence of 
nearly contemporaneous optical emission in most GRBs would be a strong argument 
against FBM. In MHD models with cr > 1 reverse shock is weak, while in EMM reverse 
shock is absent altogether. 
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Since possible observation of reverse shock emission may play a decisive role in 
validating a GRB theory, next we discuss briefly properties of the re verse sh ock ex pected 
within a framework of FBM (for more extensive discussion see, e.g. Sari & Piran. 11999 : 
Kobavashi 2000). In the framework of FBM both reverse shock and internal shocks 


which produce prompt y-ray emission originate in the same fluid and have similar 
properties (e.g. being weakly relativistic). One can naturally expect that microphysical 
properties of particle accelerating, being complicated and poorly understood, are the 
same for the same type of shocks. Thus, the conventionally introduced quantities like 
6b (magnetic held value with respect to equipartition), e e (electron energy density with 
respect to equipartition), and 7 m i n (minimum Lorentz factor of accelerated electrons) 
must be the same for both cases. As a result, for any given burst, observations of the 
prompt emission can be used to predict properties of a corresponding optical hash. 

The amount of energy dissipated in reverse shock is comparable to the energy 
dissipated in the forward shock and to the total GRB energy (iSari fe P i ran . [1999). 
The principal difference between prompt and optical emitting electrons is the radii 
of emission and ratio of cooling to expansion time scales. In a framework of FBM, 
prompt emission is generated at distances tgrb rs./ 2 T 2 0 5tc ~ 10 12 - 10 13 cm (St ~ 0.01 
s is variability time scale of the central source and, within a framework of FBM, of 
the prompt emission), while reverse shock emission is typically generated at distances 
R ~ 2 ct $ r q ~ 10 16 cm (seen at observer time t a b s ~ t s ] we concentrate on a simple 
so called ”thick shell” case). Using conventional hreball parameterization for minimum 
Lorentz factor of accelerated particles y m j n ~ e e (m p /m e )T s where is the shock Lorentz 
factor, and parameterizing energy density of magnetic held in the plasma rest frame to 
ion energy density, pc 2 = L / (Ar^Y'/r 2 c) and b = yJe B \J < &'Kpc 1 , the ratio of prompt to 
reverse shock frequencies is 


ugrb r RS 

- rsj - 

Mrs Tgrb 

so that the peak of reverse shock emission occurs at ~ 1 — 10 eV. 

An important qualitative difference between prompt and optical emitting electrons 
is that the former are in fast cooling regime, while the latter are in slow cooling regime. 
The radius beyond which optically emitting electrons enter slow cooling regime, 


A r 


RS 


r RS 

r GRB 


St 


( 11 ) 


(Trs ~ l)Lm p a T u T , . 

r cool ~ e B e e - 37rc 3 m 2 r 3 - 2 x 10 L 50 cm, (12) 

(T rs - 1 ~ 1 and r 0 = 300 was assumed) is typically smaller that Rrs- As result, 
only small fraction of energy received by an optical electron Rrs/(T or cool) 0.02 is 
emitted; the rest is lost to adiabatic expansion. For a GRB of E 1 ~ 10 51 ergs, optical 
hash would have E ~ 10 49 ergs. For a typical GRB burst with fluency ~ 10 -6 erg/cm 2 , 
and duration of ~ 100 sec, this will result in optical hash of magnitude ~ 12m. Even 
if we increase the estimate of Rrs and duration of optical hash each by an order of 
magnitude, resulting optical hash would have ~ 17m. On the other hand, brightest 
bursts may reach huences ~ 10 -4 erg/cm 2 , which, according to these estimates, can 
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produce optical flash of ~ 7m. In addition, adiabati c coo ling results in flux decay oc t 
and a clear radio signal is expected (e.g. iNakar & Piran, 2001). Thus, fireball model 
makes a predictions that all GRBs must have optical flashes in the range 12m—17m, with 
some variations of few magnitudes (both brighter and dimmer) depending on particular 
properties of each burst. 

In the Swift era, not a single GRB has shown the predicted behavior. This is despite 
the fast on-board optical telescope and a large number of ground based robotic telescopes 
(RAPTOR, ROTSEE, TAROT and others). Reverse shock emission is virtually an 
unavoidable prediction of the fireball model, so that absence of predicted reverse shocks 
emission in the Swift era argues against the fireball model. Naturally, there is a number 
of ways that through which optical flashes can be sup p ressed (e.g. c ooling of optically 
emitting electrons on photons of prompt emission Beloborodov (2002), ’’thin shell case”, 
when the r everse s hock emis sion is spread over longer times, producing weaker signal 


T-9- 


McMahon et al. 


A possible explanation of an absence of clear reverse 
shock signal is that ejecta plasma is strongly magnetized. In the case when energy 
density of magnetic field dominates the total energy dens ity (cr > 1) reverse shock 
becomes very inefficient in dissipating flow energy (Keuuel. A C orpn iti, 1984!). 


Irregular opti cal flashes (like GRB 050525a ( Kl otz et al 


2005 ) and GRB 050904 


(Bo er et all 120051) ) may be produced by other mechanisms, like gamma-ray pair 


produ ction in front of the forward shock, iBeloborodovI (J20021) : iThompson fc Mad an 
(2000), or be a low energy tail of prompt emission. 


10.2. Electromagnetic model: bright early afterglows 

Fireball and electromagnetic models m ake very different prediction for the properties 


of early afterglows (see Fig. IA21 Lyutikov, [2004 ). According to EMM, at the early 


afterglow stage the Lorentz factor and peak frequency are larger (and falling with time) 
than in the FBM (constant Lorentz factor and peak frequency). Early afterglow in the 
EMM are more energetic than in FBM, Fig. \M and can blend with the prompt phase. 


10.3. Emission radius of prompt photons and early Swift afterglows 

One of the surprising early results from Swift satellite was detection of X-ray spikes 
and breaks in light curves at intermediate t imes, much longer than burst duration but 


well before the conven tional jet break (e.g. iTagliaferrief al 


2005 


Nousek et al. 


2005 


Chincarini et al., 2005). A typical behavior includes fast-slow-fast decay with transitions 
near 100 — 1000 seconds and ~ 10 4 seconds. This presents a real challenge to GRB 
models, since if the emission is seen ’’head on”, within angle 0 < 1/T, the radii at which 
these features should be produced correspond to radii much larger than deceleration 
radius. At these times most of the energy is in the forward shock which should produce 
smooth light curve. (Late time i njection or sp e cific distribution of Lorentz facto rs are 
some possibilities discussed (e.g. Izhang et al. . 120051 Lazzati fc Begelman l . 20051) . For 


discussion in the framework of the cannonball model see Dado et al. 
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The initial fast decaying part of aft erglows was argued to be a ” sideways” prompt 


emission, coming from angles 9 > 1/T (Kumar & Panaitescu, 


2000 ; 


Barthelmv et al. 


2Q05). If this interpretation is correct, one can determine emission radii of the prompt 


emission and compare them with model predictions. (We remind that FBM predicts 
radii of emission r em 2T^c5t ~ 10 12 — 10 13 cm, while EMM predicts r em < rdec ~ 10 16 
cm, Eq. 0 . If emission is generated at r em and is coming to observer from large angles, 
9 > 1/r, its delay with respect to the start of the prompt pulse is At ~ (r em /c)6 2 / 2. 
For typical observer angle 9 ~ 0.1 and first break of a light curve at At ~ 1000 
seconds, the implied emission radius is r em ~ 6 x 10 15 cm. This is at least two orders of 
magnitude larger than is assumed in the fireball model, but is close to the assumption of 
the electromagnetic model. [To be consistent with FBM and variability on short times 
scales, the Lorentz factor of the flow should be r 0 ~ 3000, but this would imply that 
emission is strongly de-boosted, r 0 # ~ 300 1.] Interpretation of light curves breaks 

at ~ 10 3 s as been due to prompt emission seen at large angles, 9 > 1/r, is inconsistent 
with the fireball model. 


10.4- Fast variability from large radii 

If prompt emission is produced at distances ~ 10 15 — 10 16 cm, how can fast variability, 
on times scales as short as milliseconds, be achieved? One possibility, is that emission 
is beamed in the outfl ow frame, for example due to relativistic motion of ” fundamental 
emitters” ( Lv utikov fc B landford, 2003). Possible origin of relativistic motion of 
” fundamental emitters” may be the fact that in case of relativistic reconnection occurring 
in p lasm a with a 7> 1 . the_outflowing matter reaches relativistic speeds with 7 out ~ a 
( Lvutikov fc Uzdenskv . hoQ.oli . Internal synchrotron emission by such jets, or Compton 
scattering of ambient photons will then be strongly beamed in the frame of the outflow. 

Consider an outflow moving with a bulk Lorentz factor T with randomly distributed 
emitters moving with respect to the shell rest frame with a typical Lorentz factor 7 
Highly boosted emitters, moving towards an observer, have Lorentz factor 7 ~ 277 T, 
so that modest values of 7 ^ ~ 5 — 10 -C T ~ 100 — 300 suffice to produce short time 
scale variability from large distances. As the burst progresses, larger angles and more 
of internal jets producing prompt emission become visible. Most of them will be seen 
from large angles > 1/7 t hi the bulk frame, producing smooth curves. Occasionally, 
a jet at large viewing angle, 9 > 1/r, but directed towards an observer will be seen, 
producing an X-ray flare. One expects a break in the light curve at At (r em /c)9 2 / 2, 
where 9 is a viewing angle (in a structured jet model, this is the angle between the 
jet axis and direction to the observer). Afterglow should start to blend with prompt 
emission at later times. In Fig. IA4I we plot an example of a prompt light curve in 
this model (Lyutikov, in prog.). The model readily explains many unusual properties of 
early afterglows: (i) X-ray flares and light curve breaks at late times, much longer than 
conventional prompt GRB duration (extended source activity is not needed!), (ii) fast 
variability, (iii) gradual softening of the spectrum, (iv) hardernig of a spectrum during 
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X-ray flares (IBurrows et all 2005). 


10.5. Observational implications of the electromagnetic model 


In this section we give a short discussion of how the main GRB phenomena are (or may 
be) explained within a framework of EMM. 


• Jet break in afterglow GRB outflows have large opening angles, but do not 
have a jet in a proper sense. Outflows are non-isotropic so an achromatic break is 
inferred when the viewing angle is 6% ~ 1/T. 

• Structured jet. The model predicts and gives a theoretical foundation for the 
“structured jet” profile of the external shock. 

• XRF flashes. Another testable prediction of the model is that much more 
numerous X-ray flashes (XRFs) should be observed, which may be coming “from 
the sides” of the expanding shell, where the flow is less energetic and the Lorentz 
boosting is weaker. In addition, the total bolometric energy inferred for XRFs 
(from observations of afterglows before radiative losses become important) should 
be comparable to the total bolometric energy of y-ray bursts. Generally, the 
distributions of parameters of XRFs should continuously match those of GRBs. 

• Weak thermal precursor. If a fraction 1/cr ~ 0.01 — 0.1 of the magnetic energy is 
dissipated near the source, this should produce a thermal precursor with luminosity 
~ 0.01 — 0.1 of the main GRB burst. 


• Hard-soft evolution . The trend of GRB spectra to evolve from hard to soft 
during a pulse is explained as a synchrotron radiation in an expanding flow with 
magnetic field decreasing with radius B oc y/L/r (later in a pulse emission is 
produced further out where magnetic held is weaker, so that the peak energy will be 
lower; this is similar to ”radius-to-frequency mapping” in radio pulsars and AGNs). 


Amati E peak — L correlation. 

luminosity, E peak 


A correlation between peak energy and total 


2001 


2 ) follows from the assumption of a fixed 
typical emission radii and fixed minimum particle energy since B ~ \fLvt- 


\'T (Amati et al 


• Variability. Variability of the prompt emission reflects the statistical properties 
of dissipation (and not the source activity as in the FBM). Magnetic fields are 
non-linear dissipative dynamical system which often show bursty behavior with 
power law PDF. [For example, solar flares show variability on a wide range of 
temporal scales, down to minutes, which are unrelated to the time scale of 22 years 
of magnetic field generation in the tachocline.] 


Prompt and aftergl ow polarization . 


Jr’rompt ana atterg 

JCoburn & Boggs . 200,4 


Willis et a,l 


Claims of high polarization 
2005 1 if confirmed, may provide a decisive 


test of GRB models (see though R utl edg e & Fox;, 2003). The best way to produce 
polarization in the range 10 % < 11 < 60% is through synchrotron emission in large 
scale magnetic fields (Lvutikov et al 2003). (Larger polarization can only be pro- 
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duced with inverse Compton mechanism, smaller polarization can be produced by 
small scale magnetic fields.) 

Large scale field structure in the ejecta emission may also be related to polarization 
of afterglows if fields from the magnetic shell are mixed in with the shocked 
circumstcllar material. In this case, the position angle should not change through the 
afterglow while if polarization is observed both in prompt and afterglow emission the 
position angle should be the same. Also, polarization should be most independent 
of the ”jet break” moment. 


11. Conclusion 


In this contribution I outlined the underlying assumptions for the “electromagnetic 
hypothesis” for ultra-relativistic GRB outflows. The most striking implications of the 
electromagnetic hypothesis is that particle acceleration in the sources is due to direct 
dissipation of electromagnetic energy rather than shocks and that the outflows are cold, 
electromagnetically dominated flows, at least until they become strongly dissipative. 

One of the major drawback of the model is that magnetic dissipation and particle 
acceleration are very complicated processes, depending crucially on the kinetic and 
geometric properties of the plasma. This situation may be contrasted with the shock 
acceleration schemes, where a qualitatively correct result for the spectrum of accelerated 
particles, a kinetic property , can be obtained from simple macroscopic considerations 
(jump conditions). Example of the Solar corona shows that despite being complicated 
magnetic dissipation is an effective mean of particle acceleration. 

I have discussed possible observational tests of the hypothesis. In particular, 
interpretation of early afterglow features as being due to prompt emission seen at 
large angles, 6 > 1/T, allows to measure radius at which prompt emission has been 
produced. Large prompt emission radii, ~ 6 x 10 15 cm seem to be inconsistent with 
the fireball model, but close to prediction of the electromagnetic model. Internal 
relativistic motion of ” fundamental emitters” assumed within EMM may also explain 
X-ray flares during early afterglow phases (without a need for long source activity). 
An important implication of the el ectromagnetic model is that supernova explosions 


may be magnetically driven as we ll (Le blanc fc Wilsonl. Il971t iBisnovatvi-Koganl. Il971 


Wheeler et al . 1 hood I Proga et adl2003[l . 
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to be named here. In preparing this contribution I am grateful to Roger Blandford, 
Tomas Janka, Davide Lazzati, Eliud Nakar, Maurice van Putten and Stephan Rosswog 
for discussions and comments. I am also indebted to Robert Mochkovitch for shearing 
his unpublished results. 
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Appendix A. Applicability of fluid approach for blast wave 


In case of extremely high Lorentz factors of the ejecta (which require even higher 
values of a than were assumed in this paper), fluid approximation for interaction of 
magnetized ejects with ISM may break down. Consider an interface between ejecta and 
the surrounding medium in its rest frame. As a particle from the surrounding medium 
enters ejecta, it starts gyrating in magnetic held. If a fraction cr/(cr + 1) of the source 
luminosity L is in the form of magnetic held, then the turn angle in rest frame in one 
dynamical time is 




a 


2e\/nL 


(A.l) 


cr + 1 c 5 / 2 m p T 3 

In order to justify fluid approximation this should be larger than unity, which requires 

V 4 /47re 2 L\ V6 


r < 


a 


c 5 m 2 


4 x 10" 


(A.2) 


M + 1 , 

for a > 1. Thus, for any T < 4 x 10 4 an ISM particle can complete a half turn on a time 
scale short if compared with the expansion time scale. In this case, in laboratory frame 
momentum of the ejecta will be given to the particles almost instantaneously. For larger 
Lorentz factors the instantaneous hydrodynamical approximation is not applicable, but 
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Figure Al. Current flow in the electromagnetic bubble. Current flow mostly along 
the axis, on the surface of the magnetic shell, along equator and close-up at the trailing 
part of the shell. Magnetic shell is preceded by the forward shock, typically r/T 2 ahead 
of it. Non-sphericity of the shell, which is of the order ~ 1/r 2 , is enhanced. 


if particles are turned by an angle larger than ~ 1/r (lager than ~ 7r in the observer 
frame) they will still be carried with the flow. Since the rest-frame magnetic field goes 
as ~ 1 /(tr(t)), approximately linearly with time (for constant T), the rotational phase 
of a particle increases only logarithmically, / uj' B dt' oc In t. Thus, it takes a very long 
time for a particle to complete one gyration and be expelled from the ejecta. In this 
case, the ejecta will be effectively loaded with ISM particles. 

Finally, for very high Lorentz factors, 

V 4 VFL 1 / 4 


r > 


a 


o + 1 


o5/4 


8 x 10 6 


'nrir, 


(A.3) 


a particle makes a turn of less than 1/T (in the ejecta frame) on a dynamical times scale. 
In this case the ejecta just passes through ISM without much interaction and without 
slowing down. 
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Figure A2. Evolution of Lorentz factors in the fireball and electromagnetic models 
for constant external density. Relative normalization of the curves depends on the 
viewing angle (in the EMM). 
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Figure A3. X-ray afterglow in the 2-10 keV energy band for a uniform medium of 
density n = 1 cm -3 and injected power Lj SO = 10 52 erg s—1. The source is active for 
50 s, the assumed redshift is z=2.5, fraction of energy in electrons e e = 0.1, fraction 
of energy in magnetic field cb = 0.001. Blue line: electromagnetic model, red line: 
fireball model (from Mochkovitch et al., in prep.) 
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seconds 


Figure A4. Prompt emission produced by emitters moving randomly in the bulk 
frame. Emission is generated within a shell of thickness t s c = 3 x 10 12 cm moving with 
r = 100 at distance r em = T 2 t s c by randomly distributed jets with random orientation 
moving with random Lorentz factors 1 < yr < 7 T,max = 5. Each emitter is active 
for random time 0 < t' em < 0.5f s cr = t P uise,ma.x in its rest frame. Homogeneous 
jet centered on an observer with opening angle 6 = 0.1, dimensionless parameters 
Nt: / {V^T,maxS) 2 = 1.2 (probability of seeing one sub-jet ’’head-on” from angles < 1/r) 
and N(ctp U i se>rna x/2) 2 /r 2 m 6' 2 t s cT = 0.19 (efficiency of energy conversion), where N is 
total number of emitters. Intensity of emission is oc where S is Doppler factor 

and a = 0.5 is spectral index. As the burst progresses, the average Doppler factor 
<5 « t s T/t and the average flux decays as t _ ( 2+a ) = t~ 2 5 in accordance with analytical 
estimates ( Fenimore e.t al. (Lyutikov, in prog.) 






































